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Primordial fluctuations in the relative number densities of particles, or isocurvature perturbations, 
are generally well constrained by cosmic microwave background (CMB) data. A less probed mode 
is the compensated isocurvature perturbation (CIP), a fluctuation in the relative number densities 
of cold dark matter and baryons. In the curvaton model, a subdominant field during inflation later 
sets the primordial curvature fluctuation In some curvaton-decay scenarios, the baryon and cold 
dark matter isocurvature fluctuations nearly cancel, leaving a large CIP correlated with (. This 
correlation can be used to probe these CIPs more sensitively than the uncorrelated CIPs considered 
in past work, essentially by measuring the squeezed bispectrum of the CMB for triangles whose 
shortest side is limited by the sound horizon. Here, the sensitivity of existing and future CMB 
experiments to correlated CIPs is assessed, with an eye towards testing specific curvaton-decay sce¬ 
narios. The planned CMB Stage 4 experiment could detect the largest CIPs attainable in curvaton 
scenarios with more than 3 ct significance. The significance could improve if small-scale CMB po¬ 
larization foregrounds can be effectively subtracted. As a result, future CMB observations could 
discriminate between some curvaton-decay scenarios in which baryon number and dark matter are 
produced during different epochs relative to curvaton decay. Independent of the specific motivation 
for the origin of a correlated CIP perturbation, cross-correlation of CIP reconstructions with the 
primary CMB can improve the signal-to-noise ratio of a CIP detection. For fully correlated CIPs 
the improvement is a factor of ~2—3. 

PACS numbers: 98.70.Vc, 95.35.-l-d, 98.80.-k, 98.80.Cq 


I. INTRODUCTION 


The measured cosmic microwave background (CMB) 
anisotropy power spectra [ 1 , 2 ] are consistent with adia¬ 
batic primordial fluctuations, initial conditions for which 
the relative particle number densities are spatially con¬ 
stant. Adiabatic perturbations arise in the simplest in¬ 
flationary models, where a single field drives inflation 
and sets the amplitude of perturbations in all species. 
If fluctuations in particle densities or quantum numbers 
observed today are actually set by fluctuations in more 
than one field, some fraction of the primordial fluctua¬ 
tions may be isocurvature (also known as entropy) per¬ 
turbations, for which there are initial fluctuations in the 
relative particle number densities. 

The second field may be an axion [3] (and thus a dark- 
matter candidate), a curvaton [3-6] (a field that is en¬ 
ergetically subdominant during inflation but later sets 
the density fluctuations in standard-model species), or 
alternatively, inflation itself may be driven by multiple 
fields with different couplings to standard-model parti¬ 
cles [3]. Fluctuations very similar to isocurvature fluctu¬ 
ations may also arise in topological-defect models [5, 7-9]. 

Isocurvature fluctuations are defined by the entropy 
fluctuation 
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between a species i and the photons ( 7 ); where rii de¬ 
notes the background number density of the species; drii 
its spatial fluctuation; and i G { 6 , c, 2 ^, 7 }, where b de¬ 
notes baryons, c denotes cold dark matter (CDM), and 
v denotes neutrinos. These isocurvature modes individu¬ 
ally leave an imprint on the temperature and polarization 
power spectra of the CMB [5] and are highly constrained 
by current data [ 2 ]. 

There is one joint combination of isocurvature fluctu¬ 
ations that largely escapes constraints. If 

= 5 ,^ = 0 , ( 2 ) 

Pc 

then the density perturbations carried by the two isocur¬ 
vature modes cancel in this combination when both the 
baryons and CDM are nonrelativistic. This is called a 
compensated isocurvature perturbation (CIP). 

At linear order, CIPs only affect observables through 
the difference in the baryon and CDM pressure and hence 
on scales comparable to the baryonic Jeans length [10- 
14]. For the CMB, these scales are deep into the damping 
tail and the regime of secondary anisotropy dominance, 
as well as beyond the beam scale of any foreseeable CMB 
experiment. CIPs thus do not induce an observable effect 
on the CMB at linear order [10, 12, 15]. 

There are potentially observable signatures on the 
21 cm signature of neutral hydrogen at very high red- 
shifts (in absorption). Sufficiently sensitive measure¬ 
ments for a CIP detection, however, will require a fu¬ 
turistic space-based 21 cm experiment with a baseline 
that dwarfs that of ongoing/upcoming 21 cm efforts 
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like MWA/LOFAR/PAPER/SKA by an order of mag¬ 
nitude [13], 

On the other hand, since CIPs modulate the photon- 
baryon and baryon-CDM ratios, they do impact the CMB 
at higher order. By modulating these quantities in space, 
CIPs change the two-point correlations between CMB 
multipole moments in a way that allows their reconstruc¬ 
tion [16]. This fact was applied to the WMAP 9-year data 
in Ref. [17] to set upper limits on the CIP amplitude in¬ 
dependently of their origin. Similar limits follow from 
measurements of the gas fraction in massive galaxy clus¬ 
ters [18]. In these prior works, the CIP was not assumed 
to be correlated with the dominant adiabatic fluctuation. 
Here we consider an early-universe mechanism that gen¬ 
erates CIPs correlated with adiabatic fluctuations, yield¬ 
ing another detectable signature. 

Correlated isocurvature fluctuations arise naturally in 
the curvaton model, in which the curvaton, a subdom¬ 
inant field during inflation later seeds the observed pri¬ 
mordial curvature fluctuations [ 6 , 10, 19, 20]. As differ¬ 
ent species and quantum numbers may be generated by, 
before, or after curvaton decay, there are mismatches in 
their number densities which lead to isocurvature fluctu¬ 
ations including correlated CIPs. 

We assess the sensitivity of CMB anisotropy measure¬ 
ments to correlated CIPs generated in various curvaton- 
decay scenarios and find that a CMB Stage 4 [21] experi¬ 
ment could yield a detection of the largest such CIPs with 
more than 3 (t significance. The significance could im¬ 
prove to 11(7 if polarized foregrounds and systematics can 
be modeled sufficiently to make a cosmic-variance limited 
measurement out to multipoles of I = 4000. More gener¬ 
ally, we find that cross-correlation of CIP reconstructions 
with the primary CMB can improve the signal-to-noise 
ratio for detection of fully correlated CIPs by a factor of 
^2—3 depending on the specihc experiment. 

We establish that our reconstruction method [16] re¬ 
lies on a separate-universe (SU) approximation, limiting 
its use to angular scales L < 100. This has little im¬ 
pact for the signal-to-noise ratio of CIP searches using 
completed CMB experiments, but ultimately limits CIP 
reconstruction from nearly cosmic-variance limited future 
CMB polarization experiments. We correct numerical er¬ 
rors in the reconstruction noise curves of Ref. [16]; these 
errors are ultimately negligible on the scales where the 
SU approximation is valid. We also update CIP estima¬ 
tors to include sample variance from CMB R-mode po¬ 
larization, as well as covariance between CIP estimators 
based on off-diagonal correlations between different pairs 
of observables (e.g. temperature, R-mode polarization, 
and R-mode polarization). 

We begin in Sec. II by reviewing the predictions for 
the amplitude of isocurvature perturbations and their 
correlations with the adiabatic mode in nine curvaton- 
decay scenarios. In Sec. Ill, we examine the tools for 
CIP reconstruction and compute updated reconstruction 
noise spectra based on the methods from Ref. [16]. In 
Sec. IV we determine the sensitivity of future CMB ex¬ 


periments to curvaton-inspired correlated CIPs. We as¬ 
sess the improvements in signal-to-noise ratio made pos¬ 
sible by cross-correlating the CIP reconstruction with 
CMB temperature and polarization maps. We conclude 
in Sec. V. In Appendix A, we show that our reconstruc¬ 
tion methods are limited to CIP modes that are larger 
than the sound horizon at recombination. In Appendix B 
we discuss differences with the reconstruction results of 
Ref. [16]. 


II. CORRELATED CIPS IN CURVATON 
MODELS 


A. General considerations 

The curvaton cr is a light spectator scalar field during 
inflation and starts to oscillate when the Hubble scale 
H approaches the curvaton mass Wo- shortly before or 
after the inflaton ^ decays into radiation R. Once the 
curvaton starts to oscillate, it redshifts like matter and 
comes to contribute a larger and larger fraction of the 
energy density, thus generating curvature fluctuations [4, 
6 , 22-26]. 

In general both the curvaton a and inflaton (j) con¬ 
tribute to the curvature fluctuations on constant total 
density slicing C. Depending on how dominant the cur¬ 
vaton is when it decays, as quantified by 
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the relative contribution of inflaton and curvaton contri¬ 
butions to the total curvature varies, and is given by 


C. = Qi = iX-rD)C,4>+rDC,a, (4) 

where Qi is the curvature perturbation on constant den¬ 
sity Pi slicing or equivalently the energy density pertur¬ 
bation 5pi/2>{pi + Pi) on spatially flat slicing. When 
applied to particle components, C,i is also the particle 
number density perturbation on spatially flat slicing. 
Thus, i S {a, (j), b, c, v, 7 }. 

The curvaton can also generate isocurvature fluctua¬ 
tions [4, 22, 25, 27, 28] 

Sij = 3(^i — Q), (5) 

depending on how various particle numbers were gener¬ 
ated. If they were created before curvaton decay, then 
they inherit the inflaton’s fluctuations If they were 
generated by curvaton decay, they inherit the curvaton’s 
fluctuations Co-. If they were created from the thermal 
plasma after the curvaton decay, they inherit the total 
curvature perturbation In summary [25, 27], 

{ C 0 , before decay, 

c„, by decay, ( 6 ) 

C, after decay. 
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Once generated, these curvature fluctuations remain con¬ 
stant outside the horizon [23, 29, 30]. 

We are interested, in particular, in the baryon (b) and 
cold dark matter (c) isocurvature fluctuations around the 
time of recombination. We assume that lepton number 
is not related to curvaton physics, allowing us to neglect 
neutrino isocurvature perturbations [25, 27]. We thus do 
not distinguish between photons 7 and the total radia¬ 
tion. 

The remaining two isocurvature modes Sb^y and Scy 
can be reorganized into a CIP mode and a CDM isocur¬ 
vature mode, called the effective mode since it now carries 
all of the nonrelativistic matter isocurvature fluctuations, 
none of the baryon isocurvature fluctuations and only 
part of the CDM isocurvature fluctuations.^ Specifically, 
we split each curvature fluctuation as 

G=CP'’ + Cf, (7) 


where by definition the CIP mode satisfies the compensa¬ 
tion conditions and = 0 or equiv¬ 

alently 


^CIP— fb ^CIP 

( 8 ) 

and the effective mode carries the adiabatic fluctuations 
and CDM isocurvature fluctuations but no baryon isocur¬ 
vature fluctuations 


/-eff /-eff /- 
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Here the baryon fraction is 
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Pb + Pc 
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and we have assumed that the CIP mode is defined by 
compensation after both the baryons and CDM become 
nonrelativistic. 

We define the entropy perturbation carried by the two 
modes as 


sfj = 3(Cf - Cf), 


( 11 ) 


where X G {eff, CIP}. Equation (9) then implies that 
the effective mode carries only CDM isocurvature fluctu¬ 
ations, as 

= 0. (12) 

From these relations, we have 
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^ This split between CIPs and an effective mode was introduced 
in Ref. [10] with the opposite convention of an effective baryon 
isocurvature, rather than the CDM isocurvature mode. The lat¬ 
ter has since become standard (e.g. [2]). 


Together 5”®® and give an alternate representation 
of the isocurvature modes Scy and Sby The benefit of 
this representation is that because of the compensating 
baryon and CDM entropy fluctuations, the CIP mode 
corresponds to zero total isocurvature in nonrelativis¬ 
tic species [5^^^ = fbS^^^ + (1 - fb)S^^^ = 0 ], and is 
unmeasurable in linear theory, while carrying all of the 
modulation of the baryon-photon ratio (since = 0 ), 
thus inducing potentially observable changes to CMB 
anisotropy properties at second order [13, 18]. 

Now let us consider the values of S'®® and for the 
nine baryon, CDM isocurvature scenarios, obtained by 
specifying whether or not the baryon number and CDM 
are set before, by, or after curvaton decay. We use the 
notation {bx,Cy), where x,y G (before, by, after}, b de¬ 
notes baryon number, and c denotes CDM. 

Two curvaton-decay scenarios are of particular interest 
to CIPs. For the case when the baryon number is created 
by curvaton decay, (b = (a and CDM is created before, 
Cc = Ci/i) 

ceff r 

‘^c'y _ oJb '^D 

^CIP 

C C = 3(1 - ro)} {bhy, Cbefore), (14) 
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and for baryon number created before and CDM by the 
decay. 
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In these two cases, S®® can be made small by canceling 
the td and fb terms while leaving relatively large. 
The other cases are given in Table I. 

In all cases, the isocurvature modes are proportional to 
C(T ~ C<P- This implies that the cross-correlation between 
the curvature and isocurvature modes share a universal 
correlation amplitude regardless of the curvaton-decay 
scenario. For example, for the CIP mode 
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where 


(C(fe)C.(fc')) ={‘2nfS{k-k')P^^cAk), 
(c;(fc)c4fc')) =(27r)3^(fe - fc')n.c.(fc)> 

(5«P*(fc)C(fc')) =A7rf6{k - k')Psc..^{k), (17) 
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and we have used the fact that the curvaton and inflaton 
fluctuations are uncorrelated, 

Pci = i^-rDfPc,i,+rhPc.C.- ( 18 ) 

If either the curvaton or the inflaton dominates the total 
curvature, the CIP mode is fully correlated {R = ±1), 
as is the CDM-isocurvature mode S'®®. The sign of the 
correlation depends on the decay scenario. 

In fact, independently of curvaton domination or decay 
scenario, S®^ and S®® are always fully correlated and 
cannot be considered independently. Whereas individu¬ 
ally the CIP mode implies both a photon-baryon fluctu¬ 
ation and a CDM-baryon fluctuation, it can no longer 
be considered in isolation from the effective CDM isocur¬ 
vature mode. This can lead to counterintuitive results 
when considering other representations of the isocurva¬ 
ture modes and determining their joint observational ef¬ 
fects. 

For example, in the (^before, Cbefore) case, both modes 
are present and in fact set the total S^c = 0. Obviously, 
this scenario cannot be tested through a spatial modu¬ 
lation of the baryon-CDM ratio. Nonetheless, the joint 
set of modes can be described by a CIP mode which car¬ 
ries and a fully correlated CDM isocurvature mode 
where For the purposes of the tests in this 

paper, that is the more useful description, since CMB 
observables depend mainly on the photon-baryon mod¬ 
ulation of quantities like the sound speed and damping 
scale of the plasma. In this representation, both the adi¬ 
abatic and effective modes propagate in the presence of 
a CIP-modulated baryon-photon ratio. Similarly, there 
are cases [(^after? ^by) and (Rafter? ^before)] where ^ 0, 
but the effective and adiabatic modes do not see a CIP- 
induced modulation of the baryon-photon ratio. Ulti¬ 
mately (as we see below), the most interesting cases are 
those where the CIP mode is much larger than the ef¬ 
fective mode, due to observational bounds on the latter, 
making these subtleties largely irrelevant. 


B. Observational considerations for curvaton 
domination 


In forthcoming sections, we will consider the limit of 
fully correlated CIP modes (i? « ±1), a case that results 
if the curvaton completely dominates the total curvature 
fluctuation C « The inflaton fluctuations obey the 

usual relationship to tensor modes 


PC4.C<I> ^ 

Pec 16e’ 


(19) 


where r is the tensor-to-scalar ratio and e is the slow-roll 
parameter from inflation. By comparing with Eq. (18), 
we see that the curvaton contribution to the total curva¬ 
ture C is dominant over the inflaton contribution as long 
as r ^ 16e/(l — Even if gravitational waves are de¬ 

tected near the current upper limit of r ~ 0.1, curvaton 


curvature domination can still be a good approximation 
for sufficiently large e/(l — ru)^. The remaining inflaton 
contribution would then cause a small decorrelation of 
CIP modes which we ignore in the following sections. 

Under the assumption that the inflaton fluctuations 
are negligible, there are tight constraints on the CDM- 
isocurvature fraction S'®®/C that then limit the CIP am¬ 
plitude 

cCiP 

A=^ ( 20 ) 

in each of the nine scenarios. The two-sided 95% CL 
constraints from the Planck 2015 temperature and \ow-l 
polarization analysis of totally anticorrelated and cor¬ 
related isocurvature modes with no tensors combine to 
imply [2] 

^eff 

- 0.080 < -^ < 0.042 (TT-l-lowP). (21) 

The asymmetric errors reflect the fact that there is a mild 
preference for anticorrelated CDM isocurvature modes in 
the Planck TT data [31]. The standard adiabatic ACDM 
model predicts power in excess of the observations at 
low multipole moment which can be canceled by such a 
mode. This preference would strengthen if existing up¬ 
per limits to the amplitude of a primordial gravitational 
wave background are saturated in the future by a primor¬ 
dial i?-mode detection [32]. The preliminary Planck 2015 
high-Z polarization, however, disfavors the anticorrelated 
scenario and leads to the bounds [2] 

^eff 

- 0.028 < -^ < 0.036 (TT,TE,EE-flowP) (22) 
without tensors. 

Predictions for the various scenarios simplify in this 
curvaton-dominated limit. The largest CIP amplitude is 
obtained if baryon number is created by the decay and 
CDM before 



3 rp - fb 
ro 1 - /b ’ 




(Z>by, Cbefore)- 


(23) 


Note that the CDM mode can be anticorrelated in 
this case but only satisfies observational bounds if the 
decay fraction is tuned to near the baryon fraction 
I'D ~ fb- The observational bound on S^^/C, implies 
A « 3(1 —fb)/fb ~ 16.5. The converse case gives 


3 l-fb-TD 

C ro 1-fb ' 

3, (Z^before; ^by) 


(24) 


and again allows anticorrelation and can satisfy observa¬ 
tional bounds a rD is tuned to 1 — fb- 
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TABLE I. CIP and CDM isocurvature modes for the various curvaton-decay scenarios. The center two columns give the general 
case, where both curvaton and inflaton fluctuations contribute to the curvature fluctuations; whereas the right two columns 
give the curvaton-dominated fluctuation case. Italicized cases are ruled out by observational bounds. Bold-faced cases produce 
the largest CIP of the remaining ones and also predict anticorrelated effective CDM isocurvature and curvature modes. 
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Table I lists the other cases. For the ( 6 by,Cafter) and 
(Rafter, Cby) scenarlos, S'c^/C > 0 ) and to satisfy observa¬ 
tional bounds, the CIP amplitude is either proportion¬ 
ately small (|yl| ^ 10“^) or vanishing, respectively. The 
(^before; Cafter) and (l^after; ^before) scenarios cannot Sat¬ 
isfy observational bounds on S'®® /C and are hence ruled 
out. For the simultaneous scenarios, (^beforej Cbefore) can¬ 
not satisfy observational bounds, ( 6 by, Cby) predicts small 
CIP modes, and (&after, Cafter) predicts no isocurvature 
modes. 

In summary, the two cases that produce substantial 
CIP modes are the ( 6 by, Cbefore) and ( 6 before 7 Cby) scenar¬ 
ios, which predict A « 3(1 — fb)/fb ~ 16.5 and A = —3, 
respectively. Interestingly, these are also the only two 
scenarios where the CDM isocurvature mode can cancel 
the excess low multipole power in the Planck TT data. 


III. CIP RECONSTRUCTION 

CIPs leave observable imprints on the CMB. In this 
section, we review the method for CIP reconstruction 
introduced in Ref. [16] and point out the limitations 
imposed by its use of a separate-universe approxima¬ 
tion, as further illustrated in Appendix A. Reconstruc¬ 


tion methods and results are general, and do not depend 
on whether or not CIPs are generated in a curvaton sce¬ 
nario. Cross-correlations between the CIPs and the adi¬ 
abatic mode do depend on the model. Though the tech¬ 
niques again do not depend on the level of correlation, 
for results we assume fully correlated CIPs as appropriate 
for a curvaton-dominated scenario. 


In Sec. Ill A, we discuss the separate-universe response 
of off-diagonal short-wavelength CMB two-point correla¬ 
tions to the presence of a long-wavelength CIP mode. 
This response is calculated by varying background cos¬ 
mological parameters, as shown in Sec. IIIB. Each two- 
point correlation function represents a noisy measure¬ 
ment of the CIP mode which we combine to form the 
minimum variance estimator in Sec. IIIC. We discuss 
its noise properties in Sec. HID. For correlated CIP 
modes motivated by the curvaton scenarios of Sec. II, 
the reconstruction can be correlated with the CMB fields 
themselves to enhance detection, as discussed further in 
Sec. IV. 
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A. Separate-universe approximation 

The CIP mode represents a modulation of the 
baryon-photon ratio that is compensated by CDM so as 
to cancel its purely gravitational effects. Consequently, 
it leaves no imprint on the CMB to linear order. At 
second order, other modes, including the dominant adia¬ 
batic mode, propagate on a perturbed background where 
quantities such as the photon-baryon sound speed and 
damping scale are spatially modulated. A fixed CIP 
mode breaks statistical homogeneity and hence statis¬ 
tical isotropy in the CMB, and so the CIP can be re¬ 
constructed from the correlations between different CMB 
temperature and polarization multipoles that it induces. 

In Ref. [16], an approximation for characterizing these 
off-diagonal correlations was applied, based on what 
amounts to a separate-universe approximation [33] . Since 
the CIP mode does not evolve, its impact can be charac¬ 
terized by a change in cosmological parameters, so long 
as its wavelength is sufficiently large compared with the 
scale over which the modes propagate. In Appendix A, 
we show that for CMB anisotropy shortly after recom¬ 
bination, this requires the CIP mode to be larger than 
the sound horizon at that time. This amounts to the 
limit L < 100, where L is the multipole index of the 
CIP projected onto the surface of last scattering during 
recombination. Use of the expressions derived here be¬ 
yond this domain of validity will bias the associated CIP 
estimators, an issue we discuss further in Sec. IIIC and 
Appendix A. 

For modes that satisfy this approximation, we can 
treat the CIP mode shortly after recombination as a shift 
in the background 

rif, —y -|- A), 

Uc ^ Uc - UfcA, (25) 

where 

A(fi) = = D,h), (26) 

n is the direction on the sky, and is the distance to the 
CMB last-scattering surface during recombination. This 
angular field can be decomposed into multipole moments 

^ ^iMpLM(w), (27) 

LM 

SO that the restriction on the wavelength of the CIP may 
be considered as a low-pass filter where the effects are 
characterized out to L < 100. Note that this restriction 
justifies the use of a single distance in Eq. (26) rather 
than an average over the finite width of the recombina¬ 
tion era. 

In linear theory, the impact of background parameters 
on CMB power spectra 

{Xp^,Zim) = Su'Smm'C^^, (28) 


TABLE II. The response function of the various two- 

point observables in Eq. (31). 
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l + l' + L 
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EE 

//^E,dE 1 ^E,dE\ rjrL 

(C,, -bC; )Hu, 

even 

EB 

■//^E,dE 1 /^B,dB\rrL 

-bo, )H„, 

odd 

TB 


odd 

TE 

+ Cf 

even 

BB 

/ /^B ,dB 1 /^B,dB \ ttL 

(G,, -b C; )B,,, 

even 


are characterized by transfer functions 

= IJ k'^dkT^{k)Ti^{k)P^^{k) (29) 

that are given by integral solutions to the Einstein- 
Boltzmann equations of radiative transfer. Here X and Z 
are any of the CMB temperature and polarization fields 
T, E, B. Given observational bounds on the CDM isocur¬ 
vature mode, to a good approximation we can set 5®® = 0 
when evaluating the transfer functions in Eq. (29). In the 
curvaton model with = A^, the A field is correlated 
with the CMB fields through their joint dependence on C,. 
For the case where the curvaton dominates and in a flat 
cosmology, the cross power-spectrum of A and the CMB 
fields as well as the auto power-spectrum are 

described by Eq. (29) with 

Ti^{k)=Aji{kD,) (30) 

and can be numerically evaluated in GAME [15] using 
Eq. (29). 

Note that even in the presence of CIP modes, which are 
themselves statistically isotropic, two-point correlations 
are characterized by the diagonal form of Eq. (28) as long 
as (...) is understood to be the ensemble average over 
realizations of all modes. In fact, in the curvaton model, 
the ensemble average over realizations of automatically 
includes the CIP and adiabatic modes. 

Nonetheless, it is useful to artificially separate the two 
and consider the response of CMB fields to a fixed re¬ 
alization of the CIP mode. In the SU approximation, 
this fixed CIP mode is treated as simply a change in cos¬ 
mological parameters from Eq. (25) that varies across 
the sky [17]. This variation modulates the statistical 
properties of the CMB modes according to the transfer 
functions. The utility of this split is that it exposes the 
fact that there are many pairs of CMB multipoles where 
1,1' ^ L that can be used to estimate the realization of 
A lm on our sky. We denote an average over CMB modes 
with the CIP mode fixed as (.. .)cmb- This average can 
be thought of as an average over the subset of C modes 
that are smaller in wavelength than the sound horizon in 
the presence of fixed longer-wavelength ^ modes. 

The product of the source and modulation fields in real 
space leads to a convolution in harmonic space. Hence, 









7 


it connects CMB multipole moments of different Z,m in 
the same manner as a three-point function, yielding 

= Su'Smm'C^^ 

+ Y.Ar.M&m'Su;’'^ (31) 

LM 

where 


cLM _ ( 1 '\m 

V. -*-/ 


{2L + \){2l + 1){2V + l) 
47r 


/ Z L V \ 
—m M m' j ’ 


(32) 


and the response functions are given in Table II 

with 


and 


^x,dz k^^kTiX{k)^{k)P^^{k) ( 33 ) 


Kt,= 


I L r\ 

0 0 0 ) ’ 

1 L I' \ 

2 0 -2 J ’ 


which are Wigner 3j coefficients. The response for in¬ 
trinsic B modes is new to this work and may provide 
extra information on CIP modes should they be detected 
in the future. 

In the presence of a fixed long-wavelength CIP mode, 
statistical isotropy of short-wavelength CMB fields is 
therefore broken. Statistical isotropy is of course restored 
once the full ensemble average over the random realiza¬ 
tions of the CIP mode is taken. Given the correlation of 
A, C and the T, E CMB fields, a full ensemble average in¬ 
duces a three-point correlation in the CMB which corre¬ 
lates long-wavelength modes to short-wavelength power, 
i.e. a squeezed bispectrum. 

This correlation provides a way of detecting the CIP 
mode from a noisy two-point reconstruction of A as long 
as the correlation coefficient 

riXA 

^ (34) 

(jXXfjAA 



shown in Fig. 1, remains large. The sign of the correlation 
oscillates due to acoustic oscillations in temperature and 
polarization, whereas the level of correlation depends on 
the difference in projection effects between the fields. 


B. Response of CMB anisotropies to CIP modes 

Calculating the response functions in Eq. (33) requires 
varying cosmological parameters to mimic the effect of 
the CIP within the separate- universe approximation fol¬ 
lowing Ref. [16]. While Eq. (25) provides a prescrip¬ 
tion for the main effect of compensating variations of 



FIG. 1. Correlation coefficients and between 

the CIP and the CMB temperature and polarization fields, 
respectively, as a function of multipole L for A > Q. The 
sign of the correlation oscillates due to acoustic oscillations in 
temperature and polarization, whereas the level of correlation 
depends on the difference in projection effects between the 
fields. 


the background baryon and cold dark matter densities, 
there are a number of subtleties that arise from the treat¬ 
ment of the CIP mode as an angular field A(n) rather 
than a three-dimensional field that varies along the line 
of sight. CMB temperature and polarization anisotropies 
that are generated at reionization (and thus after recom¬ 
bination) break this approximation. Eor these sources 
of CMB anisotropies, the use of the SU approximation 
is limited by the horizon scale at the given time rather 
than the sound horizon shortly after recombination. 

By varying flf, and He in the transfer functions, we im¬ 
plicitly include a reionization response to the CIP that 
depends on what other parameters are held fixed. Even 
if we assume that the reionization optical depth r is held 
fixed when varying parameters in the transfer function, 
there is still an effect on the shape of the polarization 
spectra due to the implied modulation of the redshift of 
reionization To assess the impact of reionization, we 
try two other prescriptions that attempt to remove this 
sensitivity. The first case is to simply adopt a model with 
no reionization, and the second is to neglect reionization 
in evaluating Eq. (33) and then restore it using the ana¬ 
lytic damping envelope of Ref. [34] . All three prescrip¬ 
tions yield similar (at the 10%—20% level) results for the 
sensitivity of the CMB to correlated CIPs in Sec. IV. 
We conclude that that reionization only causes a small 
ambiguity for the detectability of CIPs. If they are in 
the future detected and measured precisely, then a more 
detailed prescription will be required. For simplicity, we 
adopt here the constant r prescription. 

Similarly, the CMB fields from shortly after recombi¬ 
nation are gravitationally lensed by large-scale structure 















in the foreground. Gravitational lensing of the CMB also 
produces off-diagonal two-point correlations in the pres¬ 
ence of a fixed large-scale lensing potential. Given differ¬ 
ences in the response function, it is in principle possible 
to disentangle lensing from GIF effects internally to the 
GMB. Likewise, external delensing of the CMB can help 
remove the contamination. These topics will need to be 
addressed in the future but are beyond the scope of the 
present work. They will degrade somewhat the forecasts 
for GIF detectability due to loss of degenerate modes. 
In the following sections, we treat gravitational lensing 
effects as a source of additional Gaussian noise only and 
continue to use C; to denote the unlensed CMB power 
spectrum. 

Throughout this work, we use a flat ACDM cosmology 
consistent with the 2013 Flanck results [35].^ Here Hh = 
0.049 and He = 0.268, around which we calculate the GIF 
response, and the adiabatic scalar power spectrum with 
amplitude As = 2.215 x 10“®, spectral index ris = 0.9624, 
the reionization optical depth t = 0.0925, neutrino mass 
of a single species TOj, = 0.06 eV, and Hubble constant 
h = 0.6711, which we hold fixed. We assume that the 
tensor modes (parameterized by the tensor-to-scalar ratio 
r) are negligible, and thus that there are no intrinsic B 
modes. 


C. Minimum-variance CIP estimator 

Each pair of CMB fields provides an estimate 

of a GIF mode A lm that satisfies the triangle inequality, 
\l — I'\ < L < I + I' and M = m — m!. Any single pair, 
however, is highly noisy due to the sample variance of 
the Gaussian random CMB fields and instrument noise. 
Here, we also include the change in the power spectra 


due to lensing as an additional noise source. Including 
all noise sources means that we replace Eq. (28) with 

= 5w5m^,C^^, (35) 

where 

= (36) 

with as the measurement noise power in the sky 

maps and ^s the change in the power spectrum 

due to lensing, which we treat as noise. Furthermore, 
given a XZ field pairing, the X^,^,Zim and Zi,^,Xim 
estimators have correlated noise. Likewise, the differ¬ 
ent field pairings XZ and X'Z' are also correlated. Fol¬ 
lowing the mathematically identical lensing treatment in 
Ref. [36] (see also Ref. [37]), we optimize the weight¬ 
ing of the estimators to minimize the CMB reconstruc¬ 
tion noise. Optimal weighting generalizes the familiar 
inverse-variance weighting to inverse-covariance weight¬ 
ing for covarying estimators. Differences with the results 
of Ref. [16] are discussed in Appendix B. 

It is convenient to break the inverse-covariance weight¬ 
ing into two steps. Since we want to examine each XZ 
pair individually, we first consider the 2x2 covariance of 
its multipole pairing. We can write a general estimator 
as 

^ 1 (37) 

Iml'm' 

where enforces the triangle inequality and 

are the weights to be determined. The variance of an 
unbiased estimator due to Gaussian CMB noise becomes 



— AiMp 


)cMB 






(38) 


where we have used Wick contractions and identities for the Wigner coefficients, and 

_ (2/ -I- l)(2l' + 1) 

^iv =-^ 


(39) 


Note that the covariance of the multipole permutation gives rise to the second term. We minimize the variance by 
taking a derivative of Eq. (38) with respect to W^if,, imposing the constraint of an unbiased estimator with a Lagrange 
multiplier A, obtaining 


^xx^zzwxz* ^_^y+i'+L ^xz^xz^^xz* 


'' Ll'l 


+ a 4 ’^^ = o. 


This provides the relative weights oc g^if*, where 


XZ _ 

9lu' — 


L,xz* ^xxr<zz 




-\-L ^XZ^XZ 


I'l 


c/ 




nxx nzz nxx nzz 

Lx// Lx/ Lx/ Lx// — 


(Q 


XZqXZ 


(40) 


(41) 


Specifically, we use baseline model 2.13 from Grid_limit68. 
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Note that this takes the form of the inverse 2x2 covariance weight as expected. 

The normalization comes from the requirement that the estimator be unbiased so that 


_ atXZ 7 XZcLM 

^LM L / j •) 


lml'‘ 




XZl 


qL,XZ Xz 
9lw- 


(42) 

(43) 


The normalization factor is also the variance of the estimator itself, 


\At^-ALMf) =Ni 

/ CMB 


rXZ 


(44) 


Next, we can combine the various a = XZ field pairs to find the total minimum-variance estimator by again 
inverse-covariance weighting the individual estimators 




-l\a,0 


[Nn 




al3 


(45) 


The estimator covariance-matrix Af l is at every L a rank-2 tensor over observable pairs. The indices a and j3 take 
values over labels for pairs of observables, that is, a,j3 & {TT, EE,TE, BT, BE}. Using Eq. (42), and identities of 
Wigner coefficients, we obtain an expression for the matrix elements : 

+ (-1)'+''+^ . (46) 

w 


M^zx'z' ^ j^xzjyx'z 




The total-estimator variance is again the normalization factor It is straightforward to check that Eqs. (43) and 

(44) can be recovered from Eq. (46) by restriction to a single pair {X = X', Z = Z'). 


TABLE III. Instrument noise parameters for illnstrative ex¬ 
periments [21, 38-43]: full-width half-max (FWHM) of the 
beam (in arcmin), noise for temperature measurements 5tt 
(in fiK arcmin) and, sky fraction /sky used for cosmological 
analysis. For Planck, we indicate temperature and polariza¬ 
tion noise separately as described in the text. For reconstruc¬ 
tion, we minimum variance weight the V and W bands for 
WMAP and the 143 and 217 GHz channels for Planck. The 
CVL case is full sky and has no instrument noise by definition. 


Data 

FWHM 

Noise 

/sky 

WMAP V band 

21 

434 

0.65 

WMAP W band 

13 

409 

0.65 

Planck 143 GHz 

7.1 

37, 78 

0.65 

Planck 217 GHz 

5.0 

54, 119 

0.65 

AGTPol 

1.4 

8.9 

0.097 

SPT-3G 

1.1 

2.5 

0.06 

CMB-S4 

3.0 

1.0 

0.50 

CVL 

0.0 

0.0 

1 


D. Reconstruction noise 

For reconstruction noise forecasts from the CMB 
fields XZ^ we can use the CMB noise power spec- 
ihcations (real or projected) of various experiments. We 
parameterize it as 

( 47 ) 


with Sxz the detector noise covariance assumed to 
be zero ii X ^ Z, and ^fwhm as the full-width half¬ 
max of an approximately Gaussian beam. Table III gives 
the specifications for WMAP [38], Planck [39], ACTPol 
[40, 41], SPT-3G [42, 43], and CMB Stage 4 (CMB-S4 
henceforth) [21] experiments. The WMAP and Planck 
missions have concluded, but because WMAP data have 
not yet been used to search for correlated CIPs, and 
because no CIP reconstruction yet exists from Planck 
data, we ’’predict” in those cases as well. For all but 
the Planck case, we take = 2Sgj,. For 

Planck, not all HFI bolometers have polarization sensi¬ 
tivity. To forecast temperature noise, we use all 3S + 
and diS -b 4P bolometers from the 143 and 217 GHz 
channels respectively, where S denotes an unpolarized 
spider-web bolometer and P a polarized bolometer. We 
also calculate the ideal reconstruction noise for the (zero 
instrument-noise) hypothetical cosmic-variance-limited 
(CVL) case as the ultimate limit. 

The results for CMB-S4 are shown in Fig. 2 (see Ap¬ 
pendix B for other experiments). For all cases, we gen¬ 
erate reconstruction noise curves using a maximum ob¬ 
served CMB multipole index Icmb = 2500. Beyond this 
point, foregrounds dominate the TT spectrum. To as¬ 
sess the possible effect of CMB foreground subtraction 
and lower contamination in the polarization spectra, we 
explore increasing the limit in Sec. IV. Also shown in 
Fig. 2 is the bound of ]A] < 808 on the CIP signal power 
converted from Ref. [17] for the fiducial cosmology. This 
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It is instructive to develop some intuition for the shapes 
of the reconstruction-noise curves in Fig. 2. From the ge¬ 
ometry of converting E into B alone, we can understand 
the relative slopes of between B-based and non-i?- 
based estimators of Alm- The difference is easiest to 
see in the flat-sky approximation. In order to generate 
a B mode from the CIP modulation of an E mode, the 
modulation must change the direction of the mode rela¬ 
tive to the polarization direction. In the squeezed limit 
where I, V ^ L, the flat-sky correspondence gives Fourier 
modes where I || I' and so does not generate a B mode. 
In the full-sky formulas this comes about because to good 
approximation [44] 


Htv 




sin(2^), 

cos(2(^). 


I -\- V L odd, 

I + 1' L even. 


(48) 


FIG. 2. Reconstruction noise curves for CMB-S4 ob¬ 

tained from the CIP estimator of Eq. (42) with the indi¬ 
cated pair of observables XZ. The shaded region represents 
L{L+1)C^^/(27r) signals that are excluded by WMAP 9-year 
data (Ref. [17], see Sec. IV). This bound comes from limits 
to the auto-correlation power spectrum Cf''^ of CIPs and is 
thus conservative for models with correlated CIPs. Where 
the curves intersect this bound, the estimator noise and CIP 
sample variance of the \A\ = 808 model that saturates it are 
equal. The approximate domain of validity of the SU approx¬ 
imation L < 100 is indicated by the arrow. Unless otherwise 
specified we assume, the estimators employ CMB multipoles 
up to IcMB = 2500 throughout. 

bound comes from limits to the auto-correlation power 
spectrum of CIPs. It is thus valid but could be 

further improved if A and the primordial curvature ^ are 
correlated, as we discuss further in Sec. IV. 

The CIP estimator used in all these forecasts, Eq. (42), 
is derived under the SU approximation. Naive applica¬ 
tion of Eq. (31) outside its regime of validity can thus 
significantly bias estimates of the CIP amplitude Alm- 
In the context of the simple toy model in Appendix A, we 
can compute the response exactly and estimate the re¬ 
duced sensitivity and bias beyond the SU approximation. 
Using those results, we estimate that the reconstruction 
noise curves are accurate for L < 100. Moreover, the 
sharply reduced response beyond this point in the toy 
model implies that we lose little information by simply 
restricting ourselves to these L values. 

The CIP constraints of Ref. [17] are unaffected by 
the limitations of the SU approximation, as the WMAP 
signal-to-noise ratio (S/N) for CIPs is dominated by 
scales L < 10. Planck will also not run up against the 
limitations of the SU estimator. On the other hand, 
future experiments will not be as sensitive to CIPs as 
naively estimated in the SU approach, as additional in¬ 
formation from polarization-based estimators is severely 
reduced when the L < 100 limit is imposed. We dis¬ 
cuss further the limitations to the detection of curvaton- 
generated CIPs in Sec. IV. 


where (f) is the angle between the I and I' sides of a triangle 
with side lengths {I, L, I'} upon which E also depends. In 
the squeezed limit, sm{2(p) cx L and cos{2(j)) « 1, which 
explains why oc L~^, as we can see is the 

case in Pig. 2. We also see in Fig. 2 that for L < 100, 
the noise spectra of estimators that do not involve B are 
white. This behavior is to be expected, given that the 
prominent features in the CMB carry the scale of the 
acoustic peaks so that the noise for a modulation by a 
smaller L does not depend on L. Finally, changes 
slope at roughly the sound horizon scale L « 100. This 
is because of the acoustic phase difference between T and 
E. The two terms in the response and are 

opposite in sign and roughly out of phase by half a period 
Al « 100. So for I, I' that differ by less than half a period, 
the response remains small, while for I, I' separated by 
more than a period, the response grows. Note, however, 
that the L > 100 required for the latter is the region for 
which the SU approximation fails, so this enhancement 
of the response is not of practical value. 

IV. CORRELATED CIP FORECASTS 

We now forecast constraints on the amplitude A of 
totally correlated (or anticorrelated) CIP modes for 
the various experiments in Table III. We thus restrict 
ourselves to the curvaton-dominated limit described in 
Sec. IIB. We use Fisher information matrix techniques 
in Sec. IV A, and then discuss the dependence of these 
results on various aspects of the data and assumptions in 
Sec. IV B. 


A. Fisher errors 

The Fisher information matrix Fij forecasts the 
inverse-covariance matrix of a set of parameters pi, in¬ 
cluding A, on which the auto- and cross-spectra pairs a = 
XZ oi the observed CMB and reconstructed CIP fields, 
{A, Z} e {T, E, B, A} depend. Under the assumption of 
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Gaussian statistics for these underlying fields, the Fisher 
information matrix can be approximated as 


5 ( 7 ? 




,_1\ a/3 


(D:,/3 


del 

dpj ’ 


(49) 


where Cl is the covariance matrix for an individual LM 
mode in the power spectrum estimator 


qXX',ZZ' _ qXZqX'Z' _|_ qXZ'qX'Z 


(50) 


and the /sky factor roughly accounts for the reduction 
in the 2L + 1 independent M modes due to the sky 
cut. Unless otherwise specified, we always employ the 
reconstruction-noise power spectrum for a A re¬ 

constructed from CMB multipoles up to Icmb = 2500. 



TABLE IV. 2aA detection threshold for A given all auto- and 
cross-spectra (see Sec. IV A for assumptions). Cross-spectra 
allow considerable improvement from the current bounds, and 
CMB-S4 is able to probe the largest prediction for the curva- 
ton model {A ~ 16.5) at more than 3a significance. 


Data 

2a A 

WMAP 

152 

Planck 

43.3 

ACTPol 

40.2 

SPT-3G 

38.2 

GMB-S4 

10.3 

GVL 

6.5 


In our case, we are interested in the parameter pi = A. 
Information on A is contained in the auto-spectrum of 
the reconstruction oc A^, the cross-spectra with the 
CMB T and E fields (7^^, (7|'^ oc A, and in principle 
(7£'®. We choose to neglect the information coming from 
BB, because CIP B-mode power will be swamped by the 
lensing spectrum by a factor of ^ 10^ for the WMAP9 
allowed model shown in Fig. 2 [16]. 

Since other cosmological parameters that define the 
curvature power spectrum and matter content are 
well determined, we quote 


which is the forecasted error with all other parameters 
fixed. We evaluate the derivatives in Eq. (49) at the 
fiducial ACDM model defined in Sec. IIIB. For A we 
choose the value for which a 2aA detection is possible 
which avoids problems with defining the Fisher matrix 
at A = 0 discussed below. For the multipole ranges we 
choose Lmin = max(/gj,y ,2), due to the sky cut, and 
Umax = 100, due to the breakdown of the SU approxima¬ 
tion, unless otherwise specified. 

Results for these fiducial choices are given in Table IV. 
Note that for the CMB-S4 experiment we find that the 


FIG. 3. 2a detection threshold as a function of the maxi¬ 
mum CIP multipole Lmax from all spectra for Planck, ACT- 
Pol, SPT-3G, and CMB-S4 as a function of Lmax. The shaded 
band for Lmax > 100 represents the limit of the SU approxi¬ 
mation which we take in all other results. Other parameters 
are set to the fiducial choices of Sec. IV A. 

curvaton-motivated value of A « 16.5 [realized in the sce¬ 
nario (6by, Cbefore)] Can in principle be detected at more 
than 3(7 statistical significance. Even for currently avail¬ 
able data from WMAP and Planck, the expected lim¬ 
its are substantially stronger than those determined in 
Ref. [17] for WMAP. We shall see that in large part these 
improvements are due to the addition of cross-spectra for 
correlated CIP modes. 


B. Forecast dependencies 

Now let us examine the dependence of these results 
for the CIP detection threshold on various aspects of the 
data and Fisher matrix assumptions: the maximum and 
minimum CIP multipole Umax and Uminj the maximum 
CMB multipole used in the CIP reconstruction Icmb, 
the impact of including CIP and CMB cross-correlations, 
the fiducial CIP amplitude A, and the impact of CMB 
polarization measurements on reconstruction and cross¬ 
correlation. In each case, we vary the assumptions one 
at a time from the fiducial choices in Sec. IV A. 

We begin with the multipole ranges. Recall that the 
fiducial CIP maximum multipole Umax = 100 is chosen to 
correspond to the angular scale across which sound waves 
have traveled by the end of the recombination era. As 
discussed in Appendix A, longer-wavelength CIP modes 
or smaller multipoles can be considered as SU variations 
in cosmological parameters. Since the breakdown of this 
approximation occurs within a factor of a few of this scale 
(see Fig. 7), we show in Fig. 3 the dependence of the 2a 
detection threshold on Umax- 

For the Planck experiment, whose CIP reconstruction 
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FIG. 4. 2a detection threshold as a function of the minimum 
CIP multipole Lmin from all auto- and cross-correlations for 
Planck, ACTPol, SPT-3G, and CMB-S4. Other parameters 
are set to the fiducial choices of Sec. IV A. 



FIG. 5. 2a detection threshold for the GVL case as a function 
of the maximum CMB multipole Icmb used in reconstruction. 
Other parameters are set to the fiducial choices of Sec. IV A. 


is dominated by TT estimators, the dependence on Lmax 
is very mild, consistent with the nearly white-noise spec¬ 
trum of estimators that do not involve B modes shown 
in Fig. 2. For future experiments that have good po¬ 
larization sensitivity, the noise curves of BT and BE 
can cross the others near L = 100. Were it not for the 
breakdown of the SU approximation, the implied limits 
on A would thereafter improve substantially, yielding far 
better constraints with polarization reconstruction and 
with cross-correlation than without. As discussed below, 
improvements from cross-correlation depend strongly on 
the multipole at which the signal is extracted. With the 
fiducial = 100, polarization reconstruction yields 

improvements in of order unity rather than orders 

of magnitude. 

Given that non-B-mode based CIP reconstruction has 
the highest signal-to-noise ratio at the lowest multipoles, 
it is also interesting to examine the dependence of 2a a 
on Lniin with Lniax = 100 (see Fig. 4). In addition to 
variations due to the details of the survey geometry, real 
experiments are limited by systematics, 1// noise, and 
foreground subtraction that can compromise their ability 
to probe small multipoles. As expected, the experiments 
that are the least dependent on B modes are the most 
affected by Lmin- For Planck, setting Amin = 10 degrades 
limits by a factor of 1.6, while for CMB-S4, this choice 
degrades limits by a factor of ^ 1.3. 

We have also assumed that CIP reconstruction will be 
limited to CMB multipoles smaller than Icmb = 2500 
beyond which the primary anisotropy is severely Silk 
damped. For TT-based reconstruction, foreground con¬ 
tamination will make information in higher multipoles 
difficult to extract regardless of instrument sensitivity. 
This is not necessarily the case for polarization [45], and 
so in Fig. 5, we plot the dependence on Icmb for an ideal 


CVL experiment. By Zcmb = 4000, the detection thresh¬ 
old improves to 2a a ~ 3. This level of sensitivity would 
begin testing the second largest predicted amplitude of 
the curvaton dominated scenarios A = —3 [realized in 
the scenario (^before; Cby)] and models with admixtures of 
inflaton fluctuations. Recall that these two largest cases 
are also the only ones where the effective CMB isocurva¬ 
ture mode can cancel large-angle TT power. 

Next we examine the impact of CIP cross-correlation 
with the CMB and the choice of the fiducial value of 
A. In the reconstruction noise-dominated regime, cross¬ 
correlation helps extract the signal from the noise. In 
principle, we can evaluate this improvement by compar¬ 
ing Fisher matrix errors utilizing only the auto-spectrum 
{a = AA) in Eq. (49) with the full result. There is one 
important subtlety of Fisher errors that we must address 
first. Given that oc A'^, it is clear from Eq. (49) 

that Faa oc or a a oc which diverges as A —>• 0. 
On the other hand, the corresponding limit on A^, cr ^2 
remains finite: 

(7^2 = 2|A|cta. (52) 

Of course a finite upper limit on A^ implies a finite limit 
on A as well, in spite of the Fisher estimate. The Cramer- 
Rao bound only guarantees that the Fisher estimate gives 
the best possible errors of an unbiased estimator. When 
the data only provide an upper limit, such an estima¬ 
tor can be substantially suboptimal. For this reason, we 
choose for our fiducial signal the point where A = 2a a, 
which represents a low significance detection rather than 
an upper limit. 

To ensure that our Fisher estimates of aA from auto¬ 
correlation are not misleading, consider an alternate def¬ 
inition of the error on A that corresponds to a mapping 
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of the upper limit on A^: 

a A = \/A2 + cr ^2 - \ A\ 

= v/^2 + 2|A|aA- |^|. (53) 

This quantity remains finite at >1 —>■ 0 as expected. In 
Fig. 6, we compare a a and a a of the AA auto-spectrum 
for the CMB-S4 experiment. The two estimates agree 
well as long as | A | > 2cja{A) (unshaded region). This is 
the reason we quote our primary results as the value of 
A at the detection threshold A = 2a a- 

We also show here the results obtained if only cross¬ 
spectra (TA and EA) are used, for comparison with the 
total. In the ^ 0 limit, they dominate the Fisher 

information leading a a to be finite and independent of 
the fiducial A value in this limit. 

Also interesting is the signal-dominated large-A regime 
where ^ Although the cross-spectra alone 

provide worse limits than the auto-spectrum alone, the 
total is better than what one would expect by summing 
their independent information content. This is because 
having the auto- and cross- correlation helps eliminate 
the sample variance of the Gaussian random curvature 
fluctuations Indeed, if the CIP and CMB modes were 
perfectly correlated, sample variance could be eliminated 
entirely. 

For the S4 experiment, Fig. 6 shows that the cross¬ 
spectra improve the detection threshold for A by a factor 
of 2.3. It is interesting to trace this improvement back to 
the level of correlation between the CIP and CMB modes. 
By assuming the noise-dominated regime 
as appropriate for a first detection, we can approximate 
the auto-spectra errors as 


a 


-2 

A 


AA 




(54) 


and compare this to the cross-spectrum where X € T,E: 


a 


-2 

A 


XA 


1 ^AA ^XX 

(55) 


where we recall that is the cross-correlation coeffi¬ 

cient shown in Fig. 1. As expected, cross-correlation is 
more important when the reconstruction signal-to-noise 
ratio is smaller. 

For a detection threshold A = 2a a, we can estimate 
this ratio and hence how the improvement scales with 
experimental assumptions. For the auto-spectrum this 
threshold occurs when 


^AA 


I 




(56) 


where La is a representative multipole, roughly the L 
value by which the Fisher sum accumulates half its total 
value. For this level of signal, the cross-spectra would 
give better constraints by a factor of 


o'a^Iaa 




XX 


sky ^xx 




(57) 



FIG. 6. Fisher error aA vs A for CMB-S4 from combi¬ 
nations of auto- (AA) and cross-spectra (TA, LA). Also 
shown is a a, the Fisher error implied from (j ^2 for auto¬ 
spectra from Eq. (53). The two auto-spectra analyses agree 
well for A > 2aA, which defines the regime where we can 
meaningfully compare various Fisher results. Cross-spectra 
increasingly dominate the total at low signal-to-noise ratio 
AjoA but also improve the total result at high signal-to-noise 
ratio by reducing sample variance. All other parameters are 
set according to Sec. IV A. 


Given that the correlation coefficient averaged over a suf¬ 
ficiently large range in L is always of order unity, we can 
now see that the improvement due to adding the cross¬ 
correlation depends on very few aspects of the experi¬ 
ment. For an experiment whose CIP reconstruction is 
dominated by CVL TT measurements like Planck, the 
typical CIP multipole in the signal is L < 10 and the im¬ 
provement is limited (1.9 for Planck). For an experiment 
with good polarization sensitivity and sky coverage, the 
improvement can be larger due to both the higher L out 
to which the signal can be detected and the addition of 
the EA cross-spectrum. In fact, improvements are ulti¬ 
mately limited by the SU approximation L < 100. For 
the S4 experiment, the improvement of a factor of 2.3 
from the cross-spectra comes partially from EA. With¬ 
out LA, the 2aA threshold goes from 10 to 15, and hence 
polarization plays a significant role in making a Scr de¬ 
tection of A « 16.5 possible. 

Finally, it is interesting to compare these forecasted 
constraints with existing CMB constraints to CIPs, 
which come from alone and do not apply informa¬ 
tion from cross-correlations of T/E with A. The latest 
WMAP 95% C.L. upper limit on a scale-invariant spec¬ 
trum of CIPs is L(L -b l)C'f'^ < O.OII [17]. For cor¬ 
related CIP modes, this corresponds to an upper limit 
of A^ < (808)^. The Fisher forecast from the auto¬ 
spectrum predicts that the 2a detection threshold in A^ 
is A^ = 2cr^2 = (469)^. Note that in the Fisher approxi¬ 
mation a 2cr detection in tT ^2 = 2A^ implies a a = 4A (for 
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example, the 2a threshold in A for WMAP is A = 264). 
The Fisher value in this case underestimates the actual 
errors by a factor of 3.0 in or 1.7 in A, which should 
be borne in mind when considering forecasted errors. We 
have isolated the root of this discrepancy to the differ¬ 
ence between the forecast instrument noise obtained from 
Eq. (47) and the true WMAP instrument noise at map 
level. 

The Fisher auto-result can also be compared with the 
signal-to-noise forecast in Ref. [16] for a scale-invariant 
CIP in the WMAP 7-year data release. There, to be 
maximally conservative, instrument noise was computed 
assuming a single differencing assembly. The resulting 
forecast, S/N = SOOA^j, corresponds to cr ^2 = 1.8 x 10® 
for the correlated CIP in the regime. Our 

Fisher forecast cr ^2 = 4.2 x 10^ is a factor of 4.3 lower, 
mostly due to lower values of calculated for the full 
multiple-differencing-assembly, 9-year experiment. 


V. CONCLUSIONS 

In the curvaton model, quantum fluctuations of a spec¬ 
tator field during inflation seed the primordial curva¬ 
ture perturbation after inflation and in the process can 
produce correlated isocurvature fluctuations from its de¬ 
cay. In some curvaton-decay scenarios, the usual adi¬ 
abatic and total matter isocurvature perturbations are 
accompanied by relatively unconstrained compensated 
isocurvature perturbations (CIPs) between baryons and 
dark matter. In the curvaton model, CIPs are correlated 
with the adiabatic fluctuations with amplitude given by 

= ^C. 

The most interesting (and observationally allowed) sce¬ 
narios are those where baryon number is generated by 
curvaton decay, while cold dark matter is generated be¬ 
fore; or where baryon number is generated before cur¬ 
vaton decay, and cold dark matter is generated directly 
by curvaton decay. These cases yield CIP amplitudes 
of A « 16.5 and A = —3, respectively. By modulat¬ 
ing the propagation of acoustic waves during the tightly 
coupled epoch, CIPs induce detectable off-diagonal two- 
point correlations in the CMB [16]. The correlation with 
the dominant adiabatic mode means that the cross-power 
spectrum between the CIP estimators and the CMB fields 
themselves can in principle be used to probe very small 
values of A where the auto-correlation is too noisy for 
detection. 

In this work, we obtain the expected amplitude of fully 
correlated CIPs in the different curvaton-decay scenarios 
relevant to dark matter and baryon number production. 
The sensitivity to A of seven different CMB experiments 
and the ideal cosmic-variance-limited (CVL) case was 
computed using the Fisher information and applying a 
more refined calculation of CIP reconstruction noise than 
past work [16]. We find that the validity of the implicit 
separate-universe approximation made in previous work 
[16] requires a cut L < 100 on the multipole index of the 


reconstructed CIP multipole moments Alm- While this 
cut does not affect existing limits to CIPs like Ref. [17], it 
is important for predictions of future sensitivity, partic¬ 
ularly for precise future CMB polarization experiments. 

Large-scale CMB temperature anisotropies are corre¬ 
lated with the large-scale primordial curvature pertur¬ 
bation C, and so cross-correlating the reconstructed CIP 
with CMB temperature maps can improve the detection 
threshold for A by a factor of 1.7-2.7 depending on the 
experiment. The smallest values in this range apply for 
TT-dominated experiments such as WMAP or Planck, 
and we expect that the upper limits to A from Ref. [17] 
would improve by a factor of ~ 1.7 if those CIP maps 
were cross-correlated with large-scale temperature maps. 
For a CVL experiment out to multipoles I < 2500, the 
improvement by a factor ^ 2.7 is largely independent of 
the instrument details and most sensitive to La > the mul¬ 
tipole below which the majority of the S/N comes from. 
Since the TT estimator noise has a steeper slope in L 
than the BT estimator, polarization-dominated experi¬ 
ments will naturally have S/N up to a larger La. 

The planned CMB-S4 experiment will approach the 
cosmic-variance limit for polarization. As a result, it 
could detect the A « 16.5 scenario (the largest value 
attainable in curvaton CIP scenarios) with more than 
3(7 significance. If polarized foregrounds are negligible 
or can be removed so that CIP reconstruction can be 
performed with Zcmb 4000 [45], the sensitivity to A 
of a CVL experiment will dramatically improve. In the 
cosmic-variance limit, this would allow the A « 16.5 sce¬ 
nario to be detected with ~ Her significance, and possibly 
test the second largest CIP scenario jA] = 3, as well as 
models with admixtures of inflaton fluctuations. 

A detection of fully correlated CIPs could discrimi¬ 
nate between the different curvaton-decay scenarios. The 
largest correlation A « 16.5 arises in the (5by, Cbefore) sce¬ 
nario, where the baryon number is created hy the curva¬ 
ton decay and the CDM number before the decay. This 
is the last observationally permitted scenario in which 
dark matter is produced before curvaton decay, as the 
other cases are already ruled out by the matter isocurva¬ 
ture constraints. A detection of A « 16.5 would provide 
strong support for the (5by, Cbefore) curvaton scenario, in 
which the dark matter must be produced before curva¬ 
ton decay, pointing us towards novel dark-matter pro¬ 
duction mechanisms prior to curvaton decay. This would 
also hint that baryon number generation is connected to 
the physics of a spectator field during inflation. This 
case would also predict a level of local non-Gaussianity 
of fn\ ~ 6 [46] that might be used to confirm a mea¬ 
surement from CIPs. Indeed, the Planck temperature 
(/ni = 2.5 ± 5.7, 68% C.L.) and preliminary polarization 
(/ni = 0.8 ± 5.0, 68% C.L.) constraints are already close 
to this predicted amplitude [47]. If A « 16.5 is ruled out 
by either means, we would know that in the curvaton 
model, dark matter is either directly produced by cur¬ 
vaton decay or (thermally, from the relativistic plasma) 
after curvaton decay. 
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Challenges remain for future work, in particular, a pre¬ 
cise evaluation of biases in correlated CIP measurements 
from off-diagonal correlations induced by weak gravita¬ 
tional lensing, and the generalization of the expressions 
here to models where inflaton and curvaton contributions 
to C are more comparable. Our work may pave the way 
for future CMB measurements to uncover the physics of 
curvaton decay. 


where overdots are derivatives with respect to the time 
variable rj. For the zeroth-order solution, we take Tq = 
—C/5, To = 0 as the initial condition and solve for the 
Fourier modes 

To{k, s) = -^C(fe) cos{ks), (A4) 

5 

where the sound horizon is 
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Appendix A: Wave Propagation in an 
Inhomogeneous Medium 


Given the solutions to the homogeneous equation, the 
solution for Ti is 


Ti(fe,s) 


1 d In Cg f 

5 dA J 
k'^ 

k2 _ 


§^Aik-k')Cik') 

[cos{k's) — cos(fcs)]. 


(A6) 


Note that Ti responds as an oscillator subject to an ex¬ 
ternal force given by the modulation and the unmodu¬ 
lated solution. In particular when k' = k the oscillator 
is driven at its natural frequency leading to an enhanced 
response. 

The two-point correlations with fixed modulating 
mode can likewise be expanded to first order in the mod¬ 
ulation 


{T{k)T{k'))T «(To(fc)To(fcO)T + R{k, k')AiK), (A7) 


Adiabatic acoustic waves in the GMB propagate on a 
background that is spatially modulated by the presence 
of the GIP mode Here we present a simplihed 

model of this system to test the domain of validity of 
the separate-universe (SU) approximation introduced in 
Ref. [16]. 

In this simple model, we consider an acoustic wave in 
fractional temperature fluctuations T propagating in a 
medium with temporally constant, but spatially modu¬ 
lated sound speed 


T-cf 


dine 

dA 


-A(a:) 


V^T = 0 . 


(Al) 


Note that the qualitative difference between this model 
and the SU model used in Eq. (33) is that the Taylor ap¬ 
proximation for c^(x) is employed for the medium itself, 
rather than for the observables after the acoustic mode 
has propagated through the medium. Since the medium 
is only weakly inhomogeneous, we can solve the wave 
equation by iteration. By expanding 


T = To + Ti + • • • , 


Eq. (Al) becomes 


To - cfV^To = 0, 

fi - c^V^Ti = 


din 


A(x)V^To, 


(A2) 


(A3) 


where = k + k'. The unmodulated piece is given by 


(To(k)To(k')}r = (2n)^S(k + k') cos^ks)^"^, (A8) 

25 

and the modulation response function is given by 


R{k,k')=- 


dlncl PQQ{k') 

X [cos(fc's) — cos(fcs)] -I- perm.. 


(A9) 


where the permutation refers to k k'. This response 
should be compared with the SU approximation, where 
the Taylor expansion is performed on the solution rather 
than the medium 

T^^{x, s) = To{x, s) + ^{x, s)A{x), (AlO) 


and gives off-diagonal correlations of the form of Eq. (A7) 
but with the SU response function 
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dlnCg T(;^(fc) ks 
'~dA 


dA 


perm. (All) 


sin(2fcs) -|- perm. 


Note that — cos(fcs)/5 plays the role of the transfer func¬ 
tion. 


dA 
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A comparison shows that the two are only equal in the 
limit k' —>■ k. To keep track of the differences, let us 
define 


e={k' — k)s (A12) 

and X = ks. We can rewrite the response [Eq. (A9)] as 

,A.) 

X [sin(2a: + 3e/2) — sin(e/2)] + perm. 


It is clear that the validity of the SU approximation re¬ 
quires |e| 1. It is not sufficient for the wavelength 

of the acoustic mode to be much smaller than the mod¬ 
ulating mode {k/K 3> 1). The criterion for a coherent 
driving of the oscillator is that the phase error introduced 
by e itself be small. 

Now, let us consider the relevant case for reconstruc¬ 
tion where k ^ K or x ^ e and there are many pairs 
of acoustic modes (satisfying |k| ~ |k'| ~ k) that can be 
used to measure the modulation. If is a featureless 
power law, the response can be simplified considerably, 
yielding 


R{k,k') 


d\nc^ P^dk) ks sine 


sin(2fcs -I- e). 


(AM) 


which leads to both a damping of the response and a 
decoherence in the phase for |e| > 1. 

Next, consider the impact of the reduced response on 
the estimator of the modulation mode. In the SU approx¬ 
imation, we can use Eq. (A7) to obtain the minimum- 
variance CIP estimator: 


A{K) = 


d^k T{k)T{k')R^^{k,k') 

(27r)3 



d^k [R^^{k,k')f 
(27r)3 pr^rp[^k)pTT{k') 


(A15) 


where k + k' = K and 


pTT{k) = cos^{ks)^^^ + NTT{k), (A16) 
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with Ntt as the noise power spectrum from measure¬ 
ment errors. Beyond the SU approximation, this estima¬ 
tor is biased: 


bK{Krs) = 


A{K) 


= NlP 


d^k R{k,k')R^^{k,k') 
{2ttP pTT{k)pTT{k') 


To estimate this bias, we can first determine e for each 
pair of modes that satisfies k~\-k' = K. Defining k -K = 
fikK and assuming k^ K, we have 


(A18) 


Note that e does not depend on k, but only on the angle 
of k with K. 

We can approximate the minimum variance estimator 
in Eq. (A15) by ignoring variations in the weights due 
to the k dependence of Ptt(P), in particular due to the 
unphysical zeros in power which would be filled in by 
the Doppler effect and projections in a real observable. 
Given the difference of the true- and separate-universe 
responses, this estimator would be biased as 


hK 


(f-i dfi^^sm(2x + e)[sin(2a:) -f sin(2a: -|- 2.e)])x 
{\ ill dfj,[sin{2x) + sin(2a: -I- 2e)Y)x 


2Si{2Ks) 

2Ks + sm{2Ks) ’ 


(A19) 


where denotes an average of a cycle of the oscillation 
in X. 

Near Ks « 1, the estimator becomes slightly pos¬ 
itively biased, 6^^ > 1, but quickly falls to bK < 1 
for larger values. Thus, the estimator essentially low- 
pass filters the modulation field, allowing through modes 
that are larger than the sound horizon. We show the 
behavior of bx in Fig. 7 (left), calculated using both 
Eq. (A17) (shown as a range for different noise models, 
and variety of /cmax/AT ^ 1) and the analytic approx¬ 
imation Eq. (A19). The former is nearly independent 
of assumptions and agrees very well with the latter ex¬ 
cept for a small range around it's ~ 3. Also shown is 
the result of evaluating Eq. (A17), with Ntt = 0 for 
definiteness, but with the same cycle-averaged assump¬ 
tion pTT{k) —t {pTT{k))x that was made in the ana¬ 
lytic approximation. The agreement with the analytic 
approximation shows that the overshoot around ATs ^ 3 
is due to this average. With this average, the unphysical 
zero crossings of the sound wave are eliminated in both 
Eqs. (A17) and (A19), and in that sense is closer to a 
physical model than the full calculation of Eq. (A17). In 
any case, the analytic model captures the main feature of 
the bias which is a sharp cutoff in sensitivity for Ks^ 1. 

For angular projections of the acoustic waves and the 
modulation mode, this reduced response and the corre¬ 
sponding bias is somewhat larger. In the flat-sky approx¬ 
imation, the angular modulation mode is related to the 
spatial one at an epoch s as 


A{L) = ^J = L/D,K^\)PhD^ (A20) 


where D is the distance to the observed surface, and || 
and T denote directions with respect to the line of sight. 
Defining T(Z) analogously, we again have off-diagonal cor¬ 
relations of angular moments 


{T{l)T{l')) = {To{1)To{1'))t + ^ J ^ 

xR{k,k')A{L/D,Ki\), 


e « —Ksfj,. 


(A21) 
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FIG. 7. Reconstruction bias in K (left) and in L (right) space. 
Left: Analytic expression of Eq. (A19) is in excellent agree¬ 
ment with the evaluation of Eq. (A17) using the same cycle- 
averaged replacement assumption PTT{k) —>■ {PTT(k))x. The 
gray range shows the full calculation for a wide family of Pois¬ 
son noise power spectra and fcmax. Right: Projected bias us¬ 
ing the analytic expression from Eq. (A27) with 0* = 0.01, 
showing that L < 100 is nearly unbiased. 


where L = I + 1' and 


/ dkw 

-^cos^ks) 

Pcdi/D,h) 


25112 


(A22) 


Note that here we take the angular observable as the 
value of the local temperature field on the surface. 
For the real case of CMB temperature anisotropy, the 
Doppler effect also contributes and suffers even greater 
projection effects. 

In the separate-universe approximation for R in 
Eq. (A12), Eq. (A21) reduces to 


{T{1)T{1'))t = (To(Z)To(Z'))t + (A23) 


The angular response function 
d\ncl f 

27r 4 ^ ^ 


= - 


dA 

-bperm. 


25112 


is again the derivative of the transfer function implied 
by Eq. (A22). Likewise, the SU approximation holds 
for triangles in the integrals where |e| ^ 1. Due to the 
projection, however, even if the projected 2D triangles 
in the transverse plane satisfy the analogous criteria, the 
3D triangles that compose the estimators may not, since 
K > Kj_. 

In the exact expression of Eq. (A2I), the off-diagonal 
angular multipole pairs are no longer simply proportional 
to the projected modulation mode A{L). Thus, it is 
not possible to evaluate the bias in the reconstruction 
itself. Instead, the bias appears in the auto- and cross¬ 
correlations of the reconstructed mode. Since the ex¬ 
pressions are cumbersome, we instead estimate projec¬ 
tion effects by calculating the bias that would result if 


the auto-correlation were constructed by the projection 
of biased estimators of A{K). 

For a nearly scale-invariant spectrum « 

const.). 


^AA 


£>2 J 2tt 

2n 

1yd 27r2 

27r K^P^^ 
£2 27r2 ■ 


P^^(£/£,iF||) 
" dK 


Ik^ K^yjK^-Kl 


(A25) 


Thus, the angular power spectrum gets contributions 
from modes with y = KjK^ > 1 as 


dlnCf^ , , 1 

dy ^ 2 / 2 ^^ 231 - 


(A26) 


We therefore estimate projection effects by weighting 
b'j^ accordingly: 



dyp{y)b]i{Ks 


yL9^). 


(A27) 


Here 0* = s/D \s the projected acoustic scale. In 
the ACDM cosmology, 0* « O.OI, and so £0* = 1 for 
£ = 100. The resulting projected bias is shown in Fig. 7 
(right). Clearly the bump near £0* = 1 arises from the 
bump in the three-dimensional bias plot. Fig. 7 (left). 
For this toy model of CMB acoustic waves, we see that 
£ = 100 is a scale at which to safely truncate all estima¬ 
tors of A. Since this is only a toy model, we expect that 
this estimate is only accurate to order unity corrections 
and explore sensitivity to variations in £niax in the main 
text. 


Appendix B: Improved reconstruction noise curves 

The CIP reconstruction methods introduced in 
Ref. [16] are valid as long as L < 100, as discussed in 
Appendix A. Aside from changes to precise instrumen¬ 
tal noise properties and best-fit cosmological parameters, 
the curves shown for N'['^ and in Ref. [16] are thus 
valid for all £ < 100. We find, however, that numeri¬ 
cal errors were made in Ref. [16], affecting the shape of 
at scales £ > 100; the limit to CIP reconstruction 
(£ ^ 100) imposed by the SU approximation, however, 
means that these errors are of no practical significance, 
and have no bearing on the validity of existing CMB lim¬ 
its to CIPs [17]. 

Due to an erroneous index in the code employed there, 
however, a swap took place between the indices I and I' 
when evaluating Eqs. (41) and (43) for and . 
For the last two experimental cases considered there 
[the nearly cosmic-variance-limited (CVL) EPIC mission 
concept and the actual cosmic-variance limit], the ana¬ 
lytic damping envelope reionization prescription of Ref. 
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[34], was employed, rather than the constant-r prescrip¬ 
tion employed elsewhere in Ref. [16]. For L < 100, 
we find that together these two errors lead to errors of 
< 10-2 when A e {T,E} for all exper¬ 
iments considered in Ref. [16], and so these errors are 
negligible on scales where the SU approximation is valid. 

The curves shown there for were actually nu¬ 

merically obtained using inverse-variance weighting of 
different multipole pairs, instead of the correct inverse- 
covariance weighting. Inverse-variance weights are cor¬ 
rect (and agree with inverse-covariance weights) for re¬ 
construction noise from the pairs TT, EE, TB, and EB, 
for which the observable members of a pair are either 
totally correlated or uncorrelated in the absence of an 
isotropy-breaking realization of the CIP field. For TE, 
however, the neglect of covariance between T and E leads 
to incorrect behavior. Additionally, denominators in ex¬ 
pressions for were evaluated as if Z = I' even 

when this was not the case. For L < 100, we find that 
this leads to errors of < 10“^ for all ex¬ 

periments considered in Ref. [16], and so these errors are 
negligible on scales where the SU approximation is valid. 

Calculations in Ref. [16] were sped up by using permu¬ 
tation symmetries inside sums to simplify evaluations of 
. For NI^, N'I^, and , summands were erro¬ 
neously multiplied by a factor of 1/2 when I > \L — I'\, 
but we find this leads to errors of AA'J'®/iVj® ^ 2 x 10“^ 
when L < 10 and < 5 x IQ-^ when 

10 < L < 1000, and so these errors are small on scales 
where the SU approximation is valid. 

For clarity and future reference, correct reconstruction 
noise curves NF are shown in Fig. 8 for the Planck, 
ACTPol, and SPT-3G experiments, as well as for the 
CVL case. Reconstruction noise curves for CMB-S4 are 
shown earlier in this paper, in Fig. 2. Total recon¬ 
struction noise curves are obtained using the full inverse- 
covariance weighting of different estimators [See Eqs. (45) 
and (46)] rather than the inverse-variance weighted sum 


used in Ref. [16].^ 

Using these reconstruction noise curves, we evalu¬ 
ate the signal-to-noise ratio for a detection of a scale- 
invariant spectrum of CIPs (with damping from projec¬ 
tion as in Ref. [16] for scales below the thickness of the re¬ 
combination era) using only AA auto-correlations, anal¬ 
ogously to Ref. [16], but using Planck 2013 parameter 
values [31]. Results are shown in Tables V and VI for 
a variety of experiments and for the CVL case, with 
and without the SU domain of validity imposed, and 
with/without damping from projection. We now discuss 
the results qualitatively. 

When erroneous reconstruction noise curves are used 
(with regard to the errors enumerated above), the spuri¬ 
ous improvement obtained by neglecting the SU approx¬ 
imation can be as high as ^ 50% (in the CVL case). 
If, on the other hand, the improved reconstruction noise 
curves are used, the spurious improvement falls to 5% for 
the CVL case. Put another way, if signal-to-noise ratio 
is evaluated including only modes for which the SU ap¬ 
proximation is valid, the difference between the old and 
new noise curve codes is negligible. 

If CIP projection damping is neglected (that is, a 
scale-invariant power spectrum of CIPs is assumed on all 
scales), the spurious improvement obtained by neglect¬ 
ing the SU approximation can be as high as 100% (in the 
CVL) if the erroneous reconstruction noise curves are 
used. If, on the other hand, the improved reconstruction 
noise curves are used, the spurious improvement falls to 
~ 20% for the CVL case. In all cases, the spurious im¬ 
provement in signal-to-noise (for high Umax) is caused by 
the polarization-driven flattening in the reconstruction 
noise curves at high L. 

We recommend using the reconstruction noise curves 
and tables here for future CIP forecasts in the auto- 
only case. The differences between these different sce¬ 
narios are much more dramatic if information from cross¬ 
correlations between A and T or if is included, as is the 
case in the rest of the paper. 
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TABLE V. Signal-to-noise ratio for a detection of a scale-invariant (angular) CIP power spectrum as a function of the 
variance of the baryon fraction on galaxy-cluster scales. Shown below is the constant C in the relationship SNR = Cx 1000 x 
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